Phase transitions in coupled two dimensional XY systems with spatial anisotropy 
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We study phase transitions of coupled two dimensional XY systems with spatial anisotropy and 
U(l) x I2 symmetry, motivated by spinless bosonic atoms trapped in square optical lattice on the 
metastable first excited p— level orbitals with anisotropic hopping amplitudes. The phase transitions 
of the system are generally split into an Ising transition and an XY transition, but the sequence and 
the nature of the transitions depend on the ratio between the anisotropic couplings. In the isotropic 
limit the XY variables are expected to be disordered before the Ising variables when thermal or 
quantum fluctuations are turned on gradually. In the anisotropic limit with zero perpendicular 
hoppings, the finite temperature transition is a Kosterlitz-Thouless transition driven by proliferation 
of hybrid half vortices, and the zero temperature quantum phase transition is split into a bond order 
transition and a 3D XY transition, which can be driven by the condensation of either single vortices 
or half vortices. After the condensation of half vortices the resultant state is a Mott Insulator of 
paired bosons. A small perpendicular hopping Jt, leads to a 2D Ising transition at low temperature 
and a 2+ld quantum Ising transition with a small charging energy at zero temperature. Global 
phase diagrams for both classical and quantum phase transitions are drawn. The analytical results 
obtained in this work are expected to be checked both numerically and experimentally. 

PACS numbers: 



Ever since the observation of the quantum phase tran- 
sition between the Mott Insulator and the superfluid 1 , 
the atomic systems trapped in optical lattices formed by 
laser beams have attracted a lot of attention. The po- 
tential formed by laser beams can usually be approx- 
imated by the harmonic potential. The single parti- 
cle ground state wavefunction within each well of the 
optical lattice has an approximate s— wave symmetry. 
The first excited states are three fold degenerate p— wave 
states. As discussed in reference^ and confirmed by re- 
cent experiment-, although the p— wave states are not 
the ground states, the life time for these states can be 
much longer than the average tunnelling time between 
neighboring optical lattice sites. If all of the atoms are 
pumped from the ground state to the p— wave states, they 
will establish a metastable equilibrium which can survive 
for considerable time. In previous studies, various phases 
have been identified for this system, including a super- 
fluid phase with one dimensional Z2 gauge symmetry 2,3 
on square and cubic optical lattices, a novel stripe phase 
on triangular lattice 4 , a bond algebraic liquid phased and 
a crystalline orbital order on honeycomb lattice^. Con- 
densate of p- wave paired atoms has also been considered -. 
Very recently, the long-life p— wave excited states have 
been realized experimentally, and the supcrfluid-like co- 
herence between these p— band atoms has been success- 
fully observed^. 

Orbital degeneracy also exists in materials. It is well- 
known that in transition metal oxides the d— level orbitals 
are usually partially filled, thus orbital dynamics plays 
an important role (for a recent review, see reference^) . It 
has been proposed that the quantum fluctuation of or- 
bital occupancy can lead to many interesting physics like 
"orbital liquids" 10 . However, in crystalline materials the 
quantum fluctuation which tends to disorder the orbital 



occupancy always competes with the Jahn- Teller effect 
which distorts the lattice and favors orbital ordered state 
(especially in the e g level) , therefore quantum fluctuation 
is only substantially strong in certain materials with ti g 
level degeneracy (for instance, LaTiO a 10 ' 11 ). However, in 
optical lattices formed by laser beams there is no Jahn- 
Teller effect, thus it is expected that new and interesting 
physics driven by quantum fluctuations can emerge natu- 
rally in cold atom systems on high orbital levels in optical 
lattices. 

The wave function of p x state (j) x (r) ~ x exp(—ar 2 ) 
extends further in the x than in the y and z directions, 
and so will preferentially hop to adjacent wells along the 
a;— axis. Also, p y and p z particles tend to hop in y and z 
directions respectively. If the cubic symmetry of the lat- 
tice is broken down to the planar square symmetry (for 
instance the amplitude of laser beams propagating in the 
z direction is stronger than those in the x and y direc- 
tions), the first excited p— band state is two fold degen- 
erate, with only p x and p y levels. Also, two p x particles 
can interact and convert to two p y particles, although 
single particle conversion is forbidden by the reflection 
symmetry of the square lattice. Thus the Hubbard type 
Hamiltonian describing the system can be written as 

H = V] t a b\ fiij+z - t b b\ A bi. l+y + H.c. 



H.c. 



(1) 



Here b\^ and 62,1 are annihilation operators for p x and p y 
atoms on site % respectively. The ellipses above include 
the density-density repulsive interactions generated from 
s-wave scattering between atoms. In this Hamiltonian 
(JTJ) t a > tb > 0. If the boson operators are replaced 
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by rotor operators through b Xj i ~ (— l) ia! exp(— i^) and 
b Vt i ~ (— exp(— i^a;), the Hamiltonian becomes the 
coupled anisotropic XY model 

H c = ^2 ^ J a c os{V x 9 iA ) - J b cos(Vj / ^ ij i) 

-J fc cos(V^ li2 ) - J„ cos(V v ii2 ) +7cos(26»i,i - 20 i)2 ).(2) 

In the above Hamiltonian J a ~ t a and ~ ^ are gen- 
erally unequal. Besides the classical model, we will also 
study the quantum version of this model by turning on 
the onsite repulsion energy between boson numbers n^ a , 
which is the canonical conjugate variable of 9i i0l : 

H q = H c + ^u(ni^ a - n) 2 , (3) 

a = 1, 2 denotes two orbital flavors, and n is the average 
filling of each orbital flavor on each site. 

The pair conversion term in ([1} breaks the particle 
conservation in each orbital flavor of bosons. If all the 
coefficients in ^ are finite, the symmetry of equation 
(J2J) is U(l) x Z 2 . A Z 2 variable a can be defined as 
@2 = 9i + utt/2. The U(l) phase of the system is 
X = Oi + 02- The ground state of Hamiltonian Q is a 
superfluid with both U(l) and Z 2 symmetry broken. In 
the rest of the paper we are going to study the classical 
phase transitions between the superfluid and disordered 
phase at finite temperature, as well as the quantum phase 
transition between the superfluid and the Mott Insulator 
by tuning the repulsive charging energy u. 

I. FINITE TEMPERATURE PHASE 
TRANSITIONS 

A. isotropic limit: J a = Jb 

In this section we study the finite temperature transi- 
tion in model In the isotropic limit with J a = J5, 
the system becomes the isotropic coupled XY model. 

H =J2 -^cos(V^, a )+7cos(20 M -20 2)i ). (4) 

In the absence of the pair conversion term 7, at low 
and finite temperature, the system has algebraic order of 
both exp(j#i) and exp(i# 2 ), i.e. the correlation function 
(exp(i0 Q ) ri , exp(— iO a ) r2 ) falls off as a power law with 
distance. The coupling between the two XY variables 
are always relevant in the whole algebraic phase, which 
implies that no matter how weak 7 is, at large scale 7 is 
always renormalized strong, and at long enough length 
scale one can replace 6* 2 = 9\ + cnr/2, with a = ±1. At 
this length scale, the effective Hamiltonian describing the 
coupled U(l) and Z 2 variables is 

H -ff = £ - J (! + Wi) co <°i ~ 6 J)- ( 5 ) 

<i,j> 



Now 9 = #1, and i and j are positions of blocks after 
rescaling the lattice constant. 

Many physical systems can be described by models ((4]) 
and ([5]). For instance, the J\ — J 2 antiferromagnetic XY 
model on the square lattice with nearest neighbor inter- 
action Ji and next nearest neighbor interaction J 2 breaks 
up into two square sublattices in its ground states, each 
is ordered antiferromagnetically as long as 2J 2 > J\, Ef- 
fective interaction generated by the second order spin 
wave perturbation theory favors parallel aligning between 
spins on the two sublattices, which breaks the U(l) x U (1) 
symmetry of the classical ground states down to U(l) xZ 2 
symmetr y 12 ' 13 . Another example is the fully frustrated 
XY spin model on square (or triangulr) lattice. The clas- 
sical ground states of the system break both the spin U(l) 
symmetry and the Z 2 symmetry of the chirality related 
to the rotation of spins around each unit squar e 14 ' 15 . An 
interesting fully frustrated model with J± — J 2 interac- 
tions was also studied 1 ^. 

Because in the superfluid phase both U(l) and Z 2 sym- 
metry are broken while in the high temperature disor- 
dered phase both symmetries arc restored, it is conceiv- 
able that at finite temperature there are an Ising tran- 
sition with respect to a and a Kosterlitz-Thoulcss tran- 
sition for 9 driven by the proliferation of the vortices 
of 9. A vortex of 9 is also a hybrid vortex of 9\ and 
6> 2 . A lot of efforts have been devoted to study the se- 
quence of the two transition o 17 ' 18 . Although it is still 
premature to make a conclusion, the most recent large 
size Monte-Carlo simulation suggests that there are two 
transitions at finite temperature, the Kosterlitz-Thoulcss 
transition is followed by an Ising transition at slightly 
higher temperature^, with AT C /T C ~ 0.02. Also, a the- 
oretical argument based on the fact that the cluster per- 
colation transition and the Ising transition coincide in 
the two dimensional classical Ising model indicates that 
the Ising transition should always occur at higher tem- 
perature than the Kosterlitz-Thouless transition 20 . It is 
noticed that in Hamiltonian ([5]) the XY coupling across 
an Ising domain wall is zero, thus the algebraic order 
of exp(z0) can only exist in an infinite Ising domain. Be- 
cause the density of infinite Ising domain vanishes to zero 
at the Ising critical pointSl, the Kosterlitz-Thouless tran- 
sition can only occur in the Ising ordered phase. Thus at 
the isotropic limit the Ising transition probably occurs at 
higher temperature than the XY transition. 



B. anisotropic limit: Jf, = 

In the anisotropic limit with Jb — 0, the ground state is 
superfluid order in both 9\ and # 2 with a subextensive Z 2 
degeneracy: one can perform transformation 9\ — * 9\ +7r 
(02 —> 02 + ?r) along any y (x) row, the energy of the 
state is unchanged 2 . Let us define Ising variable a as 
02 = 0i + cttt/2, and define = 0\ (although the U(l) 
phase of the system is \ = 9\ + 9%, it is more convenient 
to define 9 = 9\ to study the physics). It is noticed that 



in the anisotropic limit, without the 7 term in ([2]) both 8\ 
and 8 2 are disordered at finite temperature, therefore an 
infinitesimal 7 is not necessarily relevant. Hereafter we 
assume 7 is strong enough to guarantee that the two XY 
variables are perpendicular to each other in the tempera- 
ture range we are interested in, for instance 7 > J a ■ The 
effective classical Ising-XY coupled Hamiltonian reads 

H = ^ ~Jcos('V x 9 l ) - Ja t a l+ y cos(V y 0i). (6) 

i 

All the derivatives in the Hamiltonian are defined on the 
lattice. This model can be viewed as the XY variable 8 
coupled with an Ising gauge field background, and the 
Ising gauge field background might introduce frustration 
to the system. The lowest energy state can only be ob- 
tained if there is no frustration for 8, i.e. the Ising gauge 
flux through each plaquette is zero. The nonfrustration 
criterion generates an effective coupling between the Ising 
variables as 

Hising = -K<Ti<7 2 VzVi- (7) 

□ 

In this Hamiltonian, the summation is over all the unit 
square, and Oj, i — 1,2,3,4 are Ising variables defined 
on four corners of each unit square. This model was 
first introduced in the context of p ± ip superconductor 
array o 22,23 '. The ground state order parameter of this 
Ising Hamiltonian is the bond order bi^ — UiUi+n, this 
bond order parameter takes nonzero expectation value 
along every row and column in the whole lattice at zero 
temperature. However, any infinitesimal thermal fluctu- 
ation will disorder the bond order—. 

The finite temperature transition is driven by the £7(1) 
defects and the Z2 defects. Let us pick one ordered state 
from all the subextensive degenerate states with a = 1 
everywhere. The lowest energy U(l) defect has been 
shown to be a hybrid half vorte:^ with infinite straight 
Z2 domain walls in both x and y directions without any 
string tension, which is also shown in Fig. [T] The straight 
infinite Z2 domain walls violates the bond order (b) along 
that row and column, but as the bond order disappears at 
any finite temperature, the hybrid half vortices deconfinc 
at finite temperature. 

The lowest energy Z2 domain wall is shown in Fig. [21 
At every corner of the Z2 domain wall there is a hybrid 
half vortex attached. This implies that although the Z2 
domain wall has no string tension, the corners are inter- 
acting with each other through a logarithmic interaction. 
However, since the configuration entropy of the domain 
wall is linear with its length at finite temperature, the en- 
tropy always dominates the logarithmic interaction, thus 
the domain wall will proliferate at finite temperature. 

From Fig. Q] and Fig. [21 one can notice that across 
a Z2 domain wall along x direction, Q\ = 9 increases 
by 7r; while across a domain wall along y direction 8\ is 
unchanged on average. On the contrary, 82 gains a 7r-flip 
across a Z2 domain wall along y while no 7r-flip across Z2 
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FIG. 1: An isolated hybrid half vortex, combined with two 
infinite straight Z2 domain walls. The green (red) arrows 
stand for 61 (0 2 ). 

domain walls along x. Thus the U(l) phase \ = $1 + $2 
is disordered after the proliferation of Z2 domain walls, 
but since the superfluid stiffness is not affected by the Z2 
domain wall, 1\ is still algebraically ordered. 

There is another type of Z2 defect, which is a Z2 kink 
along one x row (see Fig. [2]). This kink defect only costs 
finite energy, and will also proliferate at finite tempera- 
ture due to thermal fluctuation (very similar to the Ising 
kinks in one dimensional classical Ising model), which 
drives the disorder of 8 along x axis. Thus one can 
conclude that at small finite temperature exp(i8) has no 
algebraic order, while exp(i2#) is algebraically ordered. 
The correlation length of 9 along x direction is similar 
to that of the one dimensional Ising model when T <C J, 
£ ~ exp(2 J/T); while along y direction the correlation 
length is zero. 

Thus in the anisotropic limit with = 0, the single fi- 
nite temperature transition is a Kosterlitz-Thouless tran- 
sition driven by the proliferation of hybrid half vortices 
shown in Fig. [TJ It is well-known that the proliferation 
of half vortices leads to a universal stiffness jump at the 
transition, which is four times as large as the stiffness 
jump driven by the proliferation of single vortices: 

A Ps /T c ~ 8/tt. (8) 

Here the stiffness is defined as the energy response to 
the same twisted boundary condition of both 8\ and 8 2 ■ 
Recently half vortex proliferation has also been predicted 
for the polar state in the two dimensional s = 1 spinor 
boson condensate 2 ^. 

The model © can be actually viewed as a "fixed 
gauge" version of the model with XY variables coupled 
with an Ising gauge field. Let us consider the following 
Hamiltonian 

Hgauge= ^ ~ J COs((j)i - . (9) 
<i,3> 
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FIG. 2: The Z2 defects. The top figure is the Z2 domain 
wall, with one hybrid half vortex attached to each corner. 
The green dashed line is the n phase mismatch of 9\ , and the 
red dashed line is the n phase mismatch of 62 ■ The Z2 domain 
walls will proliferate at infinitesimal temperature. The bot- 
tom figure is the Z2 kink which also proliferates under thermal 
fluctuation. 



Si_j = ±1 is an Ising variable defined on the link between 
sites i and j. It is well-known that this model has decon- 
fined half vortices, and the finite temperature transition 
is accompanied with a universal stiffness jump Ap s /T c ~ 
8/ iv. If now we define following gauge transforma- 

n,j,3 v =i y) j»<*a S j,j+ X > 

&i&i+y, the Hamilto- 
nian @ becomes Hamiltonian ^ . Here j x and jy are x 
and y coordinates of site j. After the gauge transforma- 
tion, an isolated half vortex with the 7r— mismatch string 
in the y direction becomes a half vortex with the string 
in the x direction. 



tions: exp(i6j) = r)j exp(i<pj), 77. 

Si,i-\-y — ^]i^li+y^i^i+y ) and frj^-j-y 



C. J a > Jb > 0, and phase diagram 

Now if a small Jh is turned up, after we introduce the 
Ising variable the Ising- XY Hamiltonian becomes 

H = - J(l + f3aiO l+3: ) cos(9i - 6 l+{i ) 

i 

-J(a i( T i+ y + (3) cos{6 l - 6 l+ y). (10) 

(3 = Jh/J a - A small string tension proportional to Jb 
is added to the Z2 domain wall in Fig. O The domain 
wall energy will compete with the configuration entropy 
of the domain wall, and the competition leads to a 2D 
Ising transition at temperature T ~ J& which separates 
an Ising ordered phase at low temperature and an Ising 
disordered phase at higher temperature. A similar Ising 
transition was found in another model with Hamiltonian 
H = E M -^cos(V^) - J a cos(2V M 6») and J b < Jj^. 
the Ising transition in this model occurs at temperature 
T ^ Jh, and the proliferation of the Ising domain walls 



implies the deconfinement of half vortices. However, in 
our case, since every corner of the Ising domain wall is 
attached with a hybrid half vortex, the deconfinement of 
half vortices would imply the proliferation of straight infi- 
nite Ising domain walls, which cost infinite energy but the 
configuration entropy is still low. Another way to look 
at this confinement is that, with small j3, the bond order 
operator — UiUi+jt always takes nonzero expectation 
value at any finite temperature, thus an infinite straight 
domain wall which destroys the bond order is not favored 
at any temperature. Thus even if the Ising domain wall 
proliferates at high temperature, the hybrid half vortices 
are still confined. The Kosterlitz-Thouless transition is 
still driven by the hybrid single vortices without any in- 
finite Z2 domain wall, which leads to a universal stiffness 
jump Ap s /T c ~ 2/tt. For the same reason, the Z2 kink 
defects are also always confined at finite temperature. 

It is noticed that after the Ising variables are disor- 
dered, the XY variable exp(z6>i) and exp(«# 2 ) cannot be 
algebraically ordered on the whole xy plane. Because of 
the symmetry between 9\ and 62, the algebraic order of 
exp(i6>i) would imply the same type of order of exp(i#2), 
which leads to the order (at least algebraic order) of Ising 
variable a ~ iexp(i8i — 1Q2). Also, since \ — @i + ^2 al- 
ways gains 7r-flip across the Ising domain wall, the U(l) 
phase x is also disordered after the Ising transition, while 
2\ is still ordered. The global finite temperature phase 
diagram is given by Fig. [3] When Jb <g; J a , the Ising 
transition occurs at temperature T ~ J^, which is lower 
than the Kosterlitz-Thouless XY transition. When Jb is 
increasing, the Ising transition finally intersects with the 
XY transition and the sequence of the transitions are re- 
versed close to the isotropic point J a = Jb- The whole 
phase diagram is symmetric for interchanging J a and Jb- 

Although the hybrid half vortices are always confined, 
as long as Jb is small enough, the experimental obser- 
vation should be similar to the anisotropic limit with 
Jb = 0. For instance, let us assume that the Kosterlitz- 
Thouless transition with Jb = occurs at temperature 
T c o, and the critical XY transition temperature with 
small Jb is T c i, then although there is no stiffness jump 
at T c o, the slope of the stiffness change is expected to be 
very steep in the window T c0 < T < T cl (Fig. 0J. The 
change of the critical temperature {T c \ — T c q)/T c o can be 
estimated as follows: With small Jb, although the hybrid 
half vortices are confined, the average distance between 
them is still large. Let us assume the average distance 
between the confined half vortices is Iq, which is much 
greater than the lattice constant. If the correlation length 
of the system is shorter than Iq, the system would behave 
as if there were deconfined hybrid half vortices. In the 
absence of Jb, the correlation length when T > T c q scales 
as £ ~ exp(y / c/t), t = T — T c q/T c q, and c is a constant 
number. The scaling formula of the correlation length 
is significantly changed with the presence of Jb only if £ 
given above is larger than l , i.e. t < c/(lni ) 2 - Thus the 
true Kosterlitz-Thouless transition driven by hybrid vor- 
tices at T c i can only occur within this small temperature 
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Jb/T (J b /u) 

FIG. 3: The phase diagram for both finite temperature clas- 
sical transitions and zero temperature quantum phase transi- 
tions. The green line is the XY transition, and the blue line 
is the Ising transition. Phase A is the ordered phase of both 
XY and Ising variables (at finite temperature there is only 
algebraic order of XY variables). Phase C is the disordered 
phase. Phase D is the phase with Ising order, and disordered 
XY variables. Phase Bi and B2 is the Ising and \ disordered 
phase, but 2\ is still ordered. 




a b 

FIG. 4: The superfluid stiffness change at the XY transition 
with zero (a) and nonzero but small (b) coupling Jt,. When 
Jb — the transition occurs at T = T c o, the stiffness jump 
is given by Ap s /T c — 8/ir; with small Jt,, although the real 
transition occurs at T c %, and the stiffness jump is Ap 3 /T c — 
2/tt, between T c o and T c \ the change of the stiffness is steep, 
thus the experimental observation is still similar to the case 
with Jt = 0. 

window. Therefore (T cl - T c0 )/T c0 ~ l/(lnZ ) 2 - 

The average distance between half vortices Iq can be 
estimated as follows: since the half vortices are attached 
to the corners of Z2 domain walls, let us consider four 
half vortices with positive vorticity at four corners of a 
rectangular Z2 domain wall. If the perimeter of the rect- 
angle is L, then the Z2 domain wall costs energy ~ JbL, 
while the repulsive logarithmic interaction between the 
half vortices gains energy ~ J a In L. Also, the number of 
configurations for rectangles is linear with L, i.e. the en- 



tropy for the rectangular Z2 domain wall is about ~ InL. 
Thus the total free energy cost for this defect configura- 
tion is F ~ JbL — (C\J a + C2T) InL, which is minimized 
when L = (Ci J a + C2T)/ Jb. Here C\ and C2 are nonuni- 
versal numbers of order 1. Thus close to the XY criti- 
cal point, the average distance between half vortices is 
^0 ~ Ja/ Jb- Combining the results in this paragraph and 
last paragraph, the change of the critical temperature 
due to Jb and the slope of the stiffness change between 
T c q and T c \ are 

T c i — T c o 1 

|g ~ [HJ a /j b )]\ (11) 

Thus when Jb / J a is small, the stiffness change still looks 
like a jump driven by half vortex proliferation. The stiff- 
ness change and jump plotted in Fig. [4] can be checked 
numerically 2 ^ 

Experimentally the ratio Jb/ J a can be tuned by chang- 
ing the depth of the wells. Superfluid stiffness jump 
was observed in superfluid helium 2 ^, and recently the 
Kosterlitz-Thouless transition has been observed in spin- 
less atomic Bose-Einstein condensates 2 ^ Our result 
about the stiffness change close to the XY transition is 
expected to be checked directly by future experiments. 



II. QUANTUM PHASE TRANSITIONS 

Now let us study the quantum phase transition at zero 
temperature by tuning the ratio between the repulsion 
energy u and the hopping energy J, and we will focus on 
the simplest case with integer average filling n, i.e. the 
system has particle-hole symmetry, more general quan- 
tum phase transitions will be discussed in another sep- 
arate work-22. The Mott insulator phase at n = 1 has 
been studied in reference^, which is a featureless Mott 
state with one p x particle and one p y particle per site. In 
the Mott insulator phase neither U(l) nor Z2 symmetry 
is broken. This quantum problem can be mapped to a 
classical problem in 3 dimensional space. In the isotropic 
case with J a = Jb = J, the corresponding three dimen- 
sional model is simply model ^ with spatial couplings in 
all three directions. In three dimensional space, although 
the argument based on the percolation of clusters^ is no 
longer applicable (because the density of infinite Ising 
domains is always finite in 3D Ising model 21 ), the U(l) 
transition is still expected to occur at lower temperature 
than the Ising transition, since the order of 9\ and 62 
implies the order of a ~ iexp[i(6i — #2)]- Thus when 
uj J is tuned from small to large, the quantum XY tran- 
sition (belonging to the 3D XY universality class) will 
occur before the quantum Ising transition (belonging to 
the 3D Ising universality class). However, the separation 
between the two transitions is expected to be very tiny 
(smaller than the two dimensional case) , because a mean 
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field Ginzburg-Landau theory for coupled XY spin sys- 
tems with the same coupling for two XY spins always 
gives rise to one single transition, and the mean field re- 
sult is expected to be more accurate in three dimensional 
systems than in two dimensional systems. 

In the anisotropic limit with Jb = 0, although the Ising 
bond order disappears with infinitesimal temperature, it 
does survive from small quantum fluctuation. Now the 
effective Hamiltonian describing the Ising variables can 
be written as 

H = J2 -Kulalalal ]T hof, (12) 

□ i 

and K ~ J , h ~ u. The bond order (bi ifl ) is nonzero 
with small h/K, because infinite order of perturbation 
of h/K is necessary in order to mix two classical ground 
states. It was conjectured that the model (|12[) has a 
quantum phase transition at J = h which separates a 
bond ordered phase and a disordered phased, and later 
on the nature of the quantum phase transition was stud- 
ied based on an exact mapping to fermionic degrees of 
freedom^. In the bond ordered phase, the hybrid half 
vortices are confined. The XY transition is driven by the 
hybrid single vortex condensation if it occurs in the bond 
ordered phase, otherwise it is driven by the condensation 
of hybrid half vortices. 

If the XY transition occurs in the bond disordered 
phase, in the dual formalism this transition can be de- 
scribed by Hamiltonian H = ^ — tcos(V A1 </) — \ai^ — 
\ a i, ai)j 4> is the phase angle of the annihilation operator 
of hybrid half vortex, and the fluxes of and a 2 are 
the boson densities ri\ and n 2 respectively. The conden- 
sate of the vortices is the Mott Insulator in the original 
boson language. One can see that after the condensation 
of hybrid half vortices, an excitation in the condensate is 
a vortex of (f>, which always involves two boson numbers 
n i + 71,2 = § ai • dl + <f 3,2 ■ dl = 2 the number of 

the vortices. Thus the condensate of hybrid half vortices 
is a Mott Insulator of paired bosons. The possible phase 
diagrams are shown in Fig. 

If a small Jf, is turned on, a quantum Ising transition 
belonging to the 3D Ising universality class is expected to 
occur at u ~ Jf,, and after the transition x = $i+#2 is dis- 
ordered while 2\ is ordered. It was shown in reference^ 
that, if an Ising system is coupled with superfluid, the 
universality class of the Ising transition is unchanged if 
the superfluid has particle hole symmetry, which is sat- 
isfied in the case we are discussing. As in the finite tem- 
perature case, with small Jb, hybrid half vortices are al- 
ways confined, and the XY transition is driven by hybrid 
single vortices condensation. After the quantum Ising 
transition, the bond order-disorder transition no longer 
exists, because the bond order parameter bi ifl — cr^a^ +fl 
always takes nonzero expectation values, hence the kink 
defect in Fig. [2] will never proliferate. The confinement 
of half vortices and Z 2 kinks can also be understood af- 
ter we map the system to a 3D classical problem, now 
a deconfined hybrid half vortex and a deconfined kink 
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FIG. 5: The phase diagram at zero temperature with increas- 
ing charging energy u along +x direction, a, J a = Jt = J, 
now the 3D XY transition occurs before the 3D Ising transi- 
tion. 6, when Jf, = 0, there are generally two transitions, one 
of which is the bond order transition at zero temperature (de- 
noted by the red closed circle). If this transition occurs first, 
the hybrid half vortices are deconfined. Phase A is the bond 
disordered phase and 2\ ordered phase. After the condensa- 
tion of hybrid half vortices, the resultant state is the Mott 
Insulator of pair bosons, c, when Jb = 0, and the bond order 
transition occurs after the XY transition, the XY transition 
is driven by the condensation of hybrid single vortices, d, 
the phase diagram with small Jf, . A quantum Ising transition 
occurs at u ~ Jb, and the bond order parameter is always 
nonzero, thus the half vortices are always confined, the XY 
transition is driven by hybrid single vortex condensation. In 
phase B after the Ising disordering, \ is disordered while 2\ 
is ordered. 



become defect lines (world lines of half vortex and kink) 
along the z direction, and the energy of these defect lines 
in the 3D space-time always dominate the entropy. For 
instance, Let us imagine a single vortex with vorticity 
+2 be split into four half vortices at (0,0,0), the four 
half vortices evolve in the 3D space time and recombine 
into a single vortex at (0, 0, z) (Fig. [6]). In any xy plane 
between (0, 0, 0) and (0, 0, z) the four half vortices are at 
four corners of a rectangle. The 3 dimensional volume 
enclosed by the world lines is an isolated Ising domain. 
The energy cost of this Ising domain is proportional to 
its boundary area, while the entropy is linear with the 
length of the half vortex world lines (for instance, the 
shape of the whole Ising domain is fully determined by 
world lines A and C in Fig. 0, i.e. the entropy is dom- 
inated by the energy and hybrid-half vortices never de- 
confinc. The confinement of Z 2 kinks can be understood 
similarly. 

Therefore we conclude that after the quantum Ising 
transition both Z 2 variable a = i exp(i#i — z# 2 ) and U(l) 
phase x = S1 + O2 are disordered, while 2% is still ordered. 
The topology of the phase diagram at zero temperature 
is the same as the finite temperature phase diagram (Fig. 
[5]). The analytical results and phase diagram Fig. [5] will 
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X 



FIG. 6: An isolated Ising domain surrounded by world lines 
of four half vortices created at (0, 0, 0) and recombine to a 
single vortex at (0, 0, z). The shape of the 3 dimensional Ising 
domain is fully determined by the world lines A and C. In 
every xy plane the for hybrid-half vortices at located at four 
corners of a rectangle. 

be checked numerically in another work2£. 

III. CONCLUSION, EXPERIMENTS AND 
OTHER PHYSICAL SYSTEMS 

The phases and phase transitions studied in this work 
can be detected experimentally by measuring the mo- 
mentum distribution function (rife) and measuring the 
momentum density correlation (nfen/-/) of the quasi- 
cquilibrium state of p— band cold atoms. In real experi- 
mental systems J& is always nonzero, let us assume Jf, is 
small and take the quantum phase transition phase dia- 
gram as an example (case d in Fig. [5]). In the phase where 
both a and are ordered, the momentum distribution at 
small momentum is 

(n k ) ~ Nk l 52 $ ~ fa °) - G) + 

G 

Nk 2 S 2 (k - (0, 7r) — (?) + 2Nk x k y S 2 (k — (n, tt) — G).(13) 

G is the basis of the reciprocal lattice. After the Ising 
transition at u ~ the sharp coherence of Q\ and 02 is 
erased, the best way to characterize the system is through 
the momentum density correlation function (rifeTifc') has 
also been proposed as a tool of measuring the many body 
effect of cold atom systems^, and this technique has 
been successfully implemented to detect crystalline den- 
sity order—. Although Q\ and 02 are not ordered after 
the Ising transitions, 20i and 202 are still ordered. Thus 
sharp coherence peaks still exist in the momentum den- 
sity correlation functions 

(n k n k ,) - (n k )(n k ,) + S 2 (k - k')(n k ) 



+c 2 NY,(k x k' x + k y k' y ) 2 S 2 (k + k' -(?) + ■■■ (14) 

G 

although one has to subtract the effect of the first two 
terms of this formula from the signal. The peaks in the 
correlation function disappear after the XY transitions. 
Thus by making use of both single particle momentum 
distribution function and momentum density correlation 
function one can detect the phase transitions studied in 
this work. 

In conclusion, we studied the phase transitions of the 
coupled anisotropic XY models. Both the classical and 
quantum phase transitions are split into an XY tran- 
sition and an Ising transition, the global phase dia- 
grams were drawn. The results obtained in this work 
can be checked numerically and hopefully by future ex- 
periments on p— band cold atom systems. Through- 
out the paper we assumed that the cubic symmetry is 
broken down to the square lattice symmetry and only 
p x and p y orbital levels are considered. If the opti- 
cal lattice has perfect cubic symmetry, then all three 
p— wave states should be taken into account. However, 
the pair conversion term between all three p— wave states 
7 {cos[2((9 1 -6 2 )]+ cos[2(0 2 - 3 )] + cos[2(0 3 -9x)}} is frus- 
trating when 7 > 0, the three terms cannot be minimized 
simultaneously, which makes the problem more compli- 
cated. 

Similar problems can be studied for atoms on even 
higher orbital levels. For instance, in cubic lattice the 
d— level state on each well is split into a three fold degen- 
erate tig level and a two fold degenerate e g level. Within 
the e g level, the d x 2_ y 2 state extends further and isotrop- 
ically in the xy plane, while the d 3z 2_ r 2 level extends in 
the z direction. Thus in the isotropic limit the prob- 
lem becomes the coupled isotropic XY model in three 
dimensional space, and as discussed in section III, the 
3D Ising transition is expected to occur at higher tem- 
perature than the 3D XY transition. On the contrary, in 
the anisotropic limit with extreme directional hoppings, 
the finite temperature model can be written as 

H = J2J2 -JWV/Ai)- Jcos(VAa) 

i fi—x : y 

+ 7 cos(26»i - 20 2 ). (15) 

If the Ising variable is introduced through 02 = 0i + (T7r/2, 
the Ising- XY model becomes 

H = J2 -Jcos(V Al i ) - Jaia i+g cos(VJi). (16) 

i /j— x,y 

Again the nonfrustration criterion generates an effective 
classical Ising Hamiltonian 

H = — J<Ji<Ji+x<Ji+z&i+x+z 

i 

-J<JiO l + yOi+~ Z Oi + y + Z. (17) 

If we define bond operator 6, = ai<r i+z , the model (fT7|) 
becomes decoupled 2D classical Ising model for new Ising 
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variable b{. Thus in the anisotropic limit there are again 
two transitions at finite temperature, one is the 3D XY 
transition, and the other is a 2D Ising transition. 

Acknowledgement 



The author thanks Joel E. Moore, Subroto Mukerjee 
and Daniel Podolsky for helpful discussions. This work is 
supported by the National Science Foundation through 
NSF grant DMR-0238760. 



1 M. Greiner and et. al., Nature 415, 39 (2002). 

2 A.Isacsson and S. M. Girvin, Phys. Rev. A 72, 053604 
(2005). 

3 W. V. Liu and C. Wu, Phys. Rev. A 74, 013607 (2006). 

4 C. Wu, W. V. Liu, J. Moore, and S. DasSarma, Phys. Rev. 
Lett 97, 190406 (2006). 

5 C. Xu and M. P. A. Fisher, Phys. Rev. B 75, 104428 (2007). 

6 C. Wu, D. Bergman, L. Balents, and S. D. Sarma, 
arXiv:cond-mat/0701788 (2007). 

7 A. B. Kuklov, Phys. Rev. Lett 97, 110405 (2006). 

8 T. Mueller, S. Foelling, A. Widera, and I. Bloch, 
arXiv:0704.2856 (2007). 

9 G. Khaliullin, Prog. Theor. Phys. Suppl 160, 155 (2005). 

10 G. Khaliullin and S. Maekawa, Phys. Rev. Lett 85, 3950 
(2000). 

11 B. Keimer and et. al., Phys. Rev. Lett 85, 3946 (2000). 

12 C. L. Henley, Phys. Rev. Lett 62, 2056 (1989). 

13 P. Chandra, P. Coleman, and A. I. Larkin, Phys. Rev. Lett 
64, 88 (1990). 

14 D. H. Lee, J. D. Joannopoulos, J. W. Negele, and D. P. 
Landau, Phys. Rev. Lett 52, 433 (1984). 

15 D. H. Lee, R. G. Caflisch, J. D. Joannopoulos, and F. Y. 
Wu, Phys. Rev. B 29, 2680 (1984). 

16 G. Franzese, V. Cataudella, S. E. Korshunov, and R. Fazio, 
Phys. Rev. B 62, R9287 (2000). 

17 J. Lee, E. Granato, and J. M. Kosterlitz, Phys. Rev. B 44, 
4819 (1991). 

18 E. Granato, J. M. Kosterlitz, J. Lee, and M. P. Nightingale, 



Phys. Rev. Lett 66, 1090 (1991). 

19 M. Hasenbusch, A. Pelissetto, and E. Vicari, J. of Stat. 
Mech p. P12002 (2005). 

20 S. E. Korshunov, Phys. Rev. Lett 88, 167007 (2002). 

21 A. Coniglio, C. R. Nappi, F. Peruggi, and L. Russo, J. of 
Phys. A 10, 205 (1977). 

22 J. E. Moore and D.-H. Lee, Phys. Rev. B 69, 104511 
(2004). 

23 C. Xu and J. E. Moore, Phys. Rev. Lett 93, 047003 (2004). 

24 S. Mukerjee, C. Xu, and J. E. Moore, Phys. Rev. Lett 97, 



25 



120406 (2006). 



D. H. Lee and G. Grinstein, Phys. Rev. Lett 55, 541 
(1985). 

26 C. Xu and S. Mukerjee, in progress (2007). 

27 D. J. Bishop and J. D. Reppy, Phys. Rev. Lett 40, 1727 
(1978). 

28 Z. Hadzibabic, P. Kruger, and et. al., Nature 441, 1118 
(2006). 

29 C. Xu, in progress (2007). 

30 H.-D. Chen, C. Fang, J. Hu, and H. Yao, Phys. Rev. B 75, 
144401 (2007). 

31 E. Frey and L. Balents, Phys. Rev. B 55, 1050 (1997). 

32 E. Altman, E. Demler, and M. D. Lukin, Phys. Rev. A 70, 
013603 (2004). 

33 S. Foiling, F. Gerbier, AWidera, and et. al., Nature 434, 
481 (2005). 



